Skyrme and Wigner crystals in graphene 
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At low-energy, the band structure of graphene can be approximated by two degenerate valleys 
(K, K') about which the electronic spectra of the valence and conduction bands have linear dispersion 
relations. An electronic state in this band spectrum is a linear superposition of states from the A 
and B sublattices of the honeycomb lattice of graphene. In a quantizing magnetic field, the band 
spectrum is split into Landau levels with level N = having zero weight on the B(A) sublattice 
for the K(K') valley. Treating the valley index as a pseudospin and assuming the real spins to 
be fully polarized, we compute the energy of Wigner and Skyrme crystals in the Hartree-Fock 
approximation. We show that Skyrme crystals have lower energy than Wigner crystals i.e. crystals 
with no pseudospin texture in some range of filling factor v around integer fillings. The collective 
mode spectrum of the valley-skyrmion crystal has three linearly-dispersing Goldstone modes in 
addition to the usual phonon mode while a Wigner crystal has only one extra Goldstone mode 
with a quadratic dispersion. We comment on how these modes should be affected by disorder and 
how, in principle, a microwave absorption experiment could distinguish between Wigner and Skyrme 
crystals. 

PACS numbers: 73.20.Qt, 73.21.-b, 73.22.Lp 



I. INTRODUCTION 

For a conventional two-dimensional electron gas 
(2DEG) created in a semiconductor heterostructure, the- 
oretical calculations show that, in the presence of a strong 
perpendicular magnetic held, a Wigner crystal (WC) 
state has lower energy then the fractional quantum Hall 
liquids for filling factors v < 1/6.5*. Transport mea- 
surements indicative of this electron crystallization has 
been reported by several groups^. These measurements 
include the observation of a strong increase in the diago- 
nal resistivity p xx , non-linear I — V characteristics, and 
broadband noise. Another series of experiments involv- 
ing microwave absorption^ have also detected a resonance 
in the real part of the longitudinal conductivity, a xx (u>), 
that has been attributed to the pinning mode of a dis- 
ordered WC. Such a resonance was observed not only 
at small filling factor v < 1/6.5 in the lowest Landau 
level but also at small filling factor in the higher Lan- 
dau levels where the formation of a quasiparticle solid is 
expected in a clean sample. In higher Landau levels, a 
study of the evolution of the pinning mode with filling 
factor reveals several transitions of the two-dimensional 
electron gas ground state from a Wigner crystal at low 
v to a series of bubble crystals with increasing number 
of electrons per lattice site as v is increased, and into 
a modulated stripe state (or anisotropic Wigner crystal) 
near half filling^. 

In a conventional 2DEG, the Landau level energy spec- 
trum is given by E n = (n + 1/2) fruj c where uj c — eB/m*c 
is the cyclotron frequency with m* the effective mass. 
Each of these levels is highly degenerate, so that a 
partially filled Landau level is dominated by electron- 
electron interactions and is expected to enter a crystal 
state. Similar physics is expected to occur in graphene. 



In a strong perpendicular magnetic field, there is also a 
series of highly degenerate Landau levels, with energies 
given by E n = sgn(n) ^^ VF yj\n\ where vf is the Fermi 



velocity. In two recent papers, Zhang and Joglekar— j£ 
have explored the possibility of Wigner crystallization in 
graphene (including bubbles and stripes) in the presence 
of a quantizing magnetic field. While the situations in the 
presence of a field are similar for the conventional 2DEG 
and graphene, in the absence of a field they are likely 
to be different. At low densities, electrons in the for- 
mer system are believed to form a Wigner crystal. How- 
ever, graphene in zero field has a gapless non-interacting 
spectrum, so that an arrangement of electrons into a 
two-dimensional lattice cannot create effective potentials 
which localize the individual electrons. This "Klein para- 
dox" physics 8 undermines the stability of the Wigner 
crystal. In a magnetic field, however, the kinetic en- 
ergy contains a series of gaps, and the physics of the 
2DEG is dominated by the Coulomb interaction alone. 
A magnetic-field induced Wigner solid is thus expected 9 . 

In the Hartree-Fock approximation (HFA), the poten- 
tial energy depends on the effective Hartree, H n (q) , and 
Fock, X n (q) interactions defined in Landau level n. To 
an excellent approximation^, these effective interactions 
for n = are identical to that of a conventional 2DEG 
so that the HFA phase diagram should be the same in 
both systems. To develop the analogy further, we use 
a pseudospin language in which the two non-equivalent 
valleys K and K' of graphene (the valence and conduc- 
tion bands in graphene touch at two inequivalent points 
K and K' = — K and four other points or valleys related 
by symmetry) are mapped to the valley pseudospin |±) 
states. If we assume that the real spins are completely 
polarized in the graphene system, then, for n = 0, the 
HF hamiltonian for the 2DEG in graphene is identical to 
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that of a conventional 2DEG with zero Zeeman coupling. 
For n 7^ 0, however, the effective interactions are not 
identical. On this basis, we can expect that pseudospin 
skyrmion crystals are possible in graphene around filling 
factor v — 1. In a conventional 2DEG, skyrmion crystals 
are restricted to the lowest Landau level only but this 
may not be the case in graphene. Indeed, recently Yang, 
Das Sarma and MacDonald have shown that skyrmions 
are the lowest-energy charged excitations in graphene for 
Landau levels up to n = 3^i. 

In this paper, we explore the possibility of Skyrmion 
crystals near integer filling in each Landau level, and 
compare their stability and collective mode properties to 
those of Wigner crystal states^. Because there is no ef- 
fective Zeeman coupling in graphene, the ground state is 
generically a crystal of merons rather than of skyrmions, 
as it is also the case in the conventional 2DEG 12 . We then 
compute the collective mode dispersions in the Wigner 
and meron crystal phases. We show that the approximate 
SU(2) symmetry of the hamiltonian leads for a Wigner 
crystal to a quadratically dispersing valley-pseudospin 
gapless mode in addition to the phonon mode present in 
the conventional 2DEG. In the meron crystal, we find in- 
stead three new gapless modes with linear dispersions in 
addition to the phonon mode. In graphene, these modes 
represent charge fluctuations between the sublattices in- 
stead of spin fluctuations as in a conventional 2DEG and 
so we speculate that they may be visible in microwave ab- 
sorption experiments. Each crystal structure, Wigner or 
Skyrme, may thus have a unique signature in microwave 
absorption, in contrast to Wigner and spin Skyrmion 
crystals in conventional 2DEG's where the absorption 
spectrum does not distinguish between these two struc- 
tures. 

This paper is organized in the following way. We re- 
view in Sees. II and III some basic properties of graphene 
and summarize our Hartree-Fock and time-dependent 
Hartree-Fock formalism for computing the phase diagram 
and collective excitations. Our numerical results are pre- 
sented in Sec. IV. We discuss in Sec. V how the collective 
modes that we find should be affected by the presence of 
disorder and speculate on their visibility in microwave 
absorption experiments. We conclude in Sec. VI. 



II. HARTREE-FOCK HAMILTONIAN 

In this section, we briefly explain the Hartree-Fock for- 
malism that we use to compute the energy and collective 
excitations of the crystal states. We start by reviewing 
the model hamiltonian for electrons in undoped graphene 
around the Fermi energy^. 

A lattice point in graphene is given by R = riiai+n2a2 
where ni,ri2 are positive or negative integers and the 
primitive vectors are chosen as ai = a/2x — V3/2y and 
a2 = ax with the two carbon atoms in the unit cell at 
positions = and r 2 — — c y where Co = a/y/3 = 
1.42 A is the separation between two adjacent carbon 



atoms. We define the carbon atoms with basis vector ri 
(r 2 ) as part of the A(B) sublattice. The tight-binding 
hamiltonian for electrons in the p z orbitals of the carbon 
atoms is given by 



H = -t. 



E 



(4h 



h.c 



(1) 



where ai (bi) is the annihilation operator for an electron 
on the A (B) sublattice of graphene at site i and the 
summation is over nearest neighbors only with hopping 
energy (between different sublattices) t « 2.8 eV. In this 
approximation, the dispersion relations for the valence 
(— ) and conduction (+) bands are given by 



E±(k) = ±t 



1 + 4 cos 2 



k x o 



(2) 



'k x a\ ( VSkyd 
+4 cos I— j cos -JL. 



1 1/2 



For undoped graphene, the Fermi level is at energy E = 
0. With our choice of orientation for the Bravais lattice, 
the positions of the two non equivalent Dirac points are at 
K = — 47r/3ax and K' = 47r/3ax. Around each of these 
points in k— space, the dispersion of the conduction and 
valence bands can be approximated by 



E (K + p) = E (K' + p) « ±hv Q |p| 



(3) 



where vo = 3coto/2H is the Fermi velocity. In the (A, B) 
basis, the hamiltonians around the Dirac points for elec- 
trons in the conduction (+) or valence (— ) band are given 
by 



H± (p) = ±Hv 



pe 



pe^ 



±i0 





(4) 



where 9 p is the angle between wavevector p and the x- 
axis. 

In the presence of a transverse magnetic field B = Bqz, 
the hamiltonian is obtained by making the Peierls substi- 
tution E (p) — > E (p + eA/hc) , where A is the vector po- 
tential of the magnetic field defined such that V x A = B. 
In terms of the covariant momentum P — — ifiV +eA/c = 
hp + eA/c, we have 



with the commutation relation 



[Px, Py] — i „ 2 



fi2 



(5) 



(6) 



where £ = y/hc/eB is the magnetic length. The original 
conical dispersions at Dirac points K and K' are now 
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split into a set of degenerate Landau levels which have 
quantized energies given by 



s y/2Hv 
E n = sgn(n) — - — v " 



(7) 



where k = 2irm/L y (with m = 0, ±1, ±2, ...) and H n (x) 
is an Hermite polynomial. 

The second quantized expression for the Coulomb in- 
teraction can be written, with the help of Eqs. (fSl fTTj) . 



where n = 0,±1,±2, ... The wavef unctions for an electron 
in these Landau levels (again in the A, B basis) are given 
by 



U 



<r|K';0,fc) 



(r|K;0,fc) 



_ -!K'T ( Vfc ( r ) 





2 ^^V(q) J dr (r|cri;ni,fci) t e lq r (r|cr 4 ;n4,fc4 



-iKr 



ho,k (r) 



(8) 
(9) 



x / dr' (r'\* 2 ; n 2 , fc 2 ) t e iqr (v'\a 3 ; n 3 , k 3 ) (13) 



for Landau level n ~ 0, and by 



1 



a/2 \sgn{ 



n)/i|„|_i, fc (r) J 
(v\K;n,k) = 1 -K-rf^WVl-U-WL) 

x 1 ' V2 V Vi.*( r ) / 



for the other levels. In the Landau gauge A = Bxy, 



h n ,k (r) = 



1/4 



-iky 



(12) 



X -H\ ('-" 2 ) 2 
xif„ ( — - — I e 



X C (Ti ,ri! ,fc! C ct 2 ,n 2 ,fe2 C<T3 !™3 ,fc3 c <r 4 ,n 4 ,fc 4 5 



where a summation over repeated indices is implied and 
where we have used er^ = ±1 for the valleys at ±K. The 
Fourier transform of the two-dimensional Coulomb inter- 
action is V (q) = 2ire 2 /q. At this point, we introduce the 

functions S^' n / (q) which we define as 



J dr (r\a; n, fc) 1 e lqr (r|c/; n' , k') (14) 
These functions S^'^, (q) are given by 



and 



+ ^2 ^ n '° e + ^''O d n D] ^"M"'! ( q ) + *n,o5 n ',0-Fo,0 (q) > 



(q) = -^2 e (M) 6 t Sffn ( n ) ( q + ct2K ) - s 9 n ("') F i»i.in'i-i (q + ct2K )] 

[<5„, e (|n'|) s 5 ti (n') F ,| B ,|_i (q + <j2K) - 5 n <, Q (|n|) s<?n (n) F H _ 1)0 (q + rfK)] 



(15) 



(16) 



where (n) is the step function and 



Fn,n'(n>ri) (q) — ( — 



V2 



cxp 



n—n I H 



(17) 



with F n! „/ ( „<„/) (q) = [-F n ', n („>„') (-q)]* • 

In our study of crystal states, we need matrix ele- 
ments of the form (G) with G ~ 2nn/ao where 
a is here the lattice constant of the Wigner or Skyrme 
crystals. We assume that the electronic density can be 
made small enough so that ao >> a, the lattice con- 
stant of graphene. Moreover, although the summations 



over G extend to infinity in the formulas below, the ex- 



ponential factor ~ e 



-G 2 £ 2 /2 



appearing in the functions 



^nn' (^) ma kes these summations rapidly convergent if 
the filling factor is not too small. We thus have an effec- 
tive cutoff value G max such that G max << K. It follows 
then that we can neglect the off diagonal matrix elements 
n' (G) * n at scatter electrons from one valley to an- 
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other since they are very small in comparison with the 
other termsi^. Essentially the same approximation was 
made in Ref. @. We also make the usual approximation 
of neglecting Landau level mixing. This approximation 
is justified since the energy of the Landau levels are given 
by Eq. (O so that the gap between the n = and n = 1 
Landau levels is thus A L . L . = V2hv F /£ « 424\/7J K 
(with B in Tesla) while the Coulomb interaction energy 
is of the order of e 2 / k£ = 130\/SK (with k = 5). It was 
recently shown numerically that Landau level mixing is 
indeed negligible^. 

With a Lande factor g ss 2, the Zeeman gap Az = 
g/J-BB — 1.34B K is however quite small in comparison 
with the Coulomb energy, and the possibility of crys- 
tal states with spin as well as valley pseudospin textures 
can also be considered. Previous studies of analogous 
bilayer 2DEG systems^ suggest that groundstates with 
real spin admixed are rather fragile with respect to Zee- 
man coupling, and it seems unlikely that such a textured 
state would be stable for this value of g. Our preliminary 
studies of the phase diagram of the combined spin and 
valley pseudospin system confirm this conclusion^. In 
what follows we assume that the electronic spin is fully 
polarized so that we need only consider the valley degree 
of freedom. 

In the Hartree-Fock approximation, our Hamiltonian 
becomes (apart from a constant term) 

H = N * E E H - W W fa)> Pn'^ (-<*) ( 18 ) 
q a, a' 

- N i> E E x « (<0 (pn a ' (q)) Pn' a (-q) . 

a, a' q 

where is the Landau level degeneracy and we have 
defined the Hartree and Fock interactions 



the average value \ p < ^ a (q)y is non zero only for q = G, 

a reciprocal lattice vector. 

The ground-state energy per electron in Landau level 
n is simply 



2i/ n 



a, a' q 



(q) 



(24) 



From Eq. (|2"Tj) . we see that the form factor So (q) for 
Landau level n = is exactly the same as for a 2DEG 
in a semiconductor quantum well or heterostructure. It 
follows that the phase diagram for graphene at low filling 
factor will be closely related to that of a conventional 
2DEG with vanishing Zeeman gapi£. 

It is very useful to map the valley degree of freedom 
into a pseudospin 1/2. Our convention is that a K state is 
pseudospin up (+) while K' is pseudospin down (— ) . In 
this language, the components of the pseudospin vector 
density P„ (q) = P n ^ x (q) x + P n>y (q) y + P n ^ z (q) z are 
given by 



P n ,x (q) 

P n,y (q) 

Pn,, (q) 



pi' 



(q) + pn' + (q) 



Pn' (q)-p n ' + (q) 
Prt -+ (q) - p«'~ (q) 



while the total density is given by 

Pn (q) = Pn" + (q) + Pn'~ (q) • 



(25) 
(26) 
(27) 

(28) 



In this language, Eq. 



becomes 



e 2 \ ( 1 



BM - UJUr" (q)5 " (_q) ' (19) 

^(q) = ^g)E^(p)^, (20) 
with the form factor 

s « (q) = 2 9 d n D (q) + F ini-i,|n|-i (q)] 

+^,0^0,0 (q) ■ (21) 
In Eq. ([IH]), we have defined the operators 

Pn a ' (q) = ^ 7 E e ^ 9a=( ' £+fe ' ) ' 2 ^^'+^in, fc c -'.«^'- 

(22) 



k.k' 



Note that, by definition, 

(pr (q = o)) = ^, 



(23) 



where v n>a — N ny(7 /N^ is the electronic filling factor of 
the n'th Landau level in the valley at crK. In a crystal, 



E 



1 



Av n \ kI ) 



e 2 \ 



^T n (q)(p„(-q))(p„(q)) (29) 



1 



^J„(q) [(P n (-q)> ■ (P n (q))] 



with 



T„(q) = 2fl„ (q) - X n (q) , 



J n (q) = (q) = 



(30) 

dye-y 2 / 2 J (yq£) , (31) 



where Jo (x) is a Bessel function. For example, the liquid 
state at v = 1 has an energy given by 



E 



ix n (0)+X n (0)|(P n (0))| 2 



(32) 



and |(P„ (0))| = 1/2. (We have taken into account the 
positive background to cancel the divergence of H n (0)). 
It follows that this liquid state is fully pseudospin polar- 
ized but the direction of polarization is arbitrary. This is 
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also true for the crystal states (see below) i.e. our hamil- 
tonian has an SU(2) symmetry. We deduce that both 
states support a pseudospin wave Goldstone mode with 
a q 2 dispersion at long wavelength. 

We remark that, in view of Eqs. ([H])-©, the pseu- 
dospin degree of freedom is equivalent to the sublattice 
degree of freedom for Landau level n = 0. This is not 
true, however, for other values of n. 



III. SINGLE AND TWO-PARTICLE GREEN'S 
FUNCTIONS 



The average values \Pn a (B)y describing the crys- 
tal states can be extracted from the Matsubara single- 
particle Green's function G n which is defined by (with 
X = k£ 2 ) 

Gr' (X, X', t) = - (Tc^x (r) cl n , x , (0)) . (33) 
Its Fourier transform is 

Gr' (G,r = 0~) = ±Y1 eri^ x ') (34) 
JV * x,x> 

xSx,x>-G v ilG%°' (X,X',t), 

so that 



pl>°(G)) = G%° (G,r = 0- 



(35) 



The equation of motion for the Green's function in the 
Matsubara formalism is obtained by using the Heisenberg 
equation 



(...) = [H HF - fiN, (. . .)] , 



(36) 



and is given by 



»(iw n -A0G£ ff (G,w„) (37) 



cr" G' 



where ui n is a fermionic Matsubara frequency and the 
matrix elements Fqq, are given by 



F^' (G, G) = H n (G - G') {Pn' a ' (G - GO) <W 

a' 

-X„(G-G0U V (G-G0). (38) 



To find the order parameters yp^ :<T {G)j , we solve 
the Hartree-Fock equation of motion numerically by an 
iterative method. The procedure is described in detail in 
Ref. [13. 

In order to compute the collective excitations, we de- 
fine the 16 response functions 

X ^ d (q,q';r) = -N+ (Tp a / (q,r) p c n d (-q', 0)) 

+^{p^(q)>(^ d (-qO>, (39) 

where a,b,c,d are valley indices. For a crystal, 
q = k + G and q = k + G' , where k is a vector in 
the first Brillouin zone. In the generalized random- 
phase approximation (GRPA), the equation of motion 
for x a,b,c,d ^ q /. ig given by 



E 

q" 



X(q",q';ifi»)=-D(q,q0 



(40) 

where / is the 4iV q x iV q unit matrix with 7V q the number 
of vectors G kept in the calculation and we have defined 
the matrices: 

x+,+.+,+ x+,+,-,+ x+,+,+ - x+,+,-,- 

X+.-+.+ X+.-- + X+,-,+ - x+--- 
X-,+,+,+ X-,+,-,+ X-+.+ - X- + -- 
X- ,-,+,+ X-,-,-,+ X--.+ - X-,-,-,- 



and 



(41) 



^(q.qO = ^(q,qO ( 42 ) 
- £ D (q, q'O (H (q", qO - X (q", q'^S,) 



q" 



with 



D = 



2ia(p+'+) -{p-- + )l (Pi'')!* 

(Pn'">7*-<^' 

(/V + >7* -{Pt'-)l 





7 -(/V + > 
2io (p n > 



7 



(44) 



6 



and 



/ -2iaH[(p+'+) + (p-'-)] 
+2iaX (p+.+) 



-2iaH[(p+<+) + (p-<-)} 







(45) 



In these equations, we adopt the conventions 
that (p a ' b ) = (p Q ' 6 (q-q')) , H (p a - b ) stands for 

^™(q-q')(p Q ' b (q-q')) I 7 = e-' iqxq '^ 2 / 2 , and a = 
sin ^ qxq ^ ^ • The Hartree and Fock interaction matri- 
ces are given by 



H 



JT„(q) ff n (q) 





H n (q) H n {q) 



(46) 



and 



X 



X n (q) 

X n (q) 

X n (q) 






Jq,q'- 



(47) 



Once the Hartree-Fock densities (p a ' b ) are calculated for 
the crystal state considered, the response functions can 
be computed using Eq. ([3D]). The collective excitations 
appear as poles of these response functions. To derive the 
dispersion relations, we follow the poles in the response 
functions as the wave vector k is varied within the first 
Brillouin zone. We consider only the low-energy modes 
in the present work. The response functions also have 
higher-energy modes corresponding to more localized ex- 
citations. 

The various response functions in Eq. (|41[) can be com- 
bined in an obvious way to give the pseudospin response 
functions x Pn , Pn , Xs x ,s x , Xs y ,s y , and xs z ,s,- 



IV. PHASE DIAGRAM AND COLLECTIVE 
MODES 

In the absence of Zeeman coupling, each Landau 
level has a fourfold degeneracy (the valley degeneracy 
combined with the usual spin doublet). For undoped 
graphene, the n — Landau level multiplet is half-filled. 
We use the notation v n £ [0,4] for the filling factor of 
each Landau level multiplet. The total filling factor is 
thus given by v — An — 2 + v n . In this work, we assume a 
finite Zeeman coupling but neglect any mixing of Landau 
level with different spins so that the phase diagram for 



u n £ [0,2] is identical with that for v n £ [2,4] . Without 
lost of generality, we consider v n £ [0, 2] from now on. 

Our procedure for solving Eq. (j37|) does not allow us 
to find the absolute ground state of the 2DEG for a given 
filling factor. Instead, we have to be content with com- 
paring the energy of different phases and finding the low- 
est one amongst them. In this study, we focus on the 
crystal states and more specifically on the Skyrme crys- 
tals. The filling factor v n is that of the partially filled 
Landau level and all filled levels below n are assumed 
inert. This procedure is valid when Landau level mixing 
is small, provided one considers only intra-Landau level 
excitations^^. We consider the following states in each 
Landau level n : 

1. Electron bubble crystal (eBCn). A triangular lat- 
tice with N c electrons per unit cell and filling factor 
v n < 0.5. More precisely, bubbles are maximal den- 
sity droplets^. 

2. Hole bubble crystal (hBCn). A triangular lattice 
with N c holes per unit cell and filling factor 0.5 < 
u n < 1. The lattice constant ao of such a crystal 
is determined by the relation 27rn^^ 2 = \v n — 1| 
where rih = N c /ea,Q is the hole density with N c the 
number of holes in the bubbles and e = V3/2 for a 
triangular lattice. Note that we find both the hBCn 
state and the meron crystal considered below to be 
lower in energy that the eBCn groundstate assumed 
in Ref. 6 in the same range of filling factors. 

3. Meron crystal (MC). A square lattice with four 
merons of charge — e/2 (if v n < 1) or e/2 (if v n > 1) 
per unit cell, equally spaced and assembled in a 
checkerboard configuration. The z component of 
the pseudospin and the vorticities alternate from 
one site to the next and two of the merons in the 
unit cell have a global pseudospin phase in the x — y 
plane which is opposite to the two others. In semi- 
conductor 2DEG's, this configuration is found for 
the (spin) skyrmion crystal when the Zeeman en- 
ergy is zeroi 2 -. This crystal is represented in Fig. 
Ill 

4. Meron pair crystal (MPC). A triangular lattice 
with four merons per unit cell. At each lattice 
site, two merons with the same value of P, and 
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vorticities but opposite values of the global phase 
are coupled together so that the pseudospins rotate 
by 4tt on a path encircling the two meron pairs. 
This configuration is represented in Fig. [5] The 
merons are not equally spaced. The possibility for 
skyrmions of o ppo site phases to form pairs was con- 
sidered in Ref. UM The MC and MPC phases arc in 
competition with each other and their energies are 
very close. We remark that, in the absence of an 
equivalent Zceman coupling, Az, a single skyrmion 
should have a size comparable to that of the sample 
size. In a lattice, this causes a strong interaction 
between skyrmions that leads to a lattice of merons 
even when the skyrmion filling factor \u n — 1| — * 0. 
It is important to notice that the limits Az — > 
and \v n — 1| — > do not commute^. If we were 
to choose \v n — 1| — > first and then Az — > 0, we 
would find instead a Skyrme crystal^. 

For the numerical calculations, we consider a filling 
factor v n G [0.1,0.9]. For v n < 0.1 or \v n - 1| < 0.1, 
the number of reciprocal lattice vectors needed in the 
calculation becomes very large and we do not get good 
convergence. This is due to the fact that the size in 
real space of the quasiparticles (electrons for v n < 0.1 or 
holes or skyrmions for \v n — 1| < 0.1) decreases so that 
more wavevectors are needed to describe them. Also, the 
hamiltonian has electron-hole symmetry around v n = 1 
so that the sequence of phase transitions found for v n > 1 
is the mirror image of that for v n < 1 with particles 
replaced by anti-particles. For example, the couterpart 
of the phase eBCl at v n = 0.2 is a hBCl at v n = 1.8 
with a filling of holes given by z/„,ft = 2.0 — 1.8 = 0.2. 
Similarly, the counterpart of a crystal of merons (with 
charge — e/2) at v n — 0.8 (with a filling of merons given 
by v n ,m — 1-0 — 0.8 = 0.2) is a crystal of anti-merons 
(charge +e/2) at v n — 1.2 with a filling of anti-merons 
given by v n>am = 1.2 - 1.0 = 0.2. 

We show in Fig. [3] the energies of different phases of 
the 2DEG in graphene for Landau level n = 0. We find 
the following sequence: eBCl for vq 6 [0.1,0.5], hBCl 
for vq E [0.5,0.55], MC for u € [0.55,0.65], and MCP 
for vq e [0.65,0.9]. As noted above, this sequence of 
transitions is the same as that calculated for a 2DEG 
in GaAs-AlGaAs quantum wells in the absence of Zee- 
man coupling because the effective interactions H (q) 
and Xq (q) are the same in both cases. To determine 
this sequence, we not only find the state with the lowest 
energy but also compute the collective mode spectrum 
in order to check that the crystal is stable. For the MC 
and MPC where the difference in energy is close to our 
numerical accuracy, the stability criteria allows us to find 
the correct ground state. 

The eBCl phase is fully pseudospin polarized but 
its energy is independent of the orientation of the 
pseudospins. The crystal thus has a full SU(2) sym- 
metry for the pseudospin. The dispersion relations 
of the two Goldstone modes of this crystal is given 
in Fig. 2J The dispersion is plotted along the path 
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FIG. 1: Pseudospin texture in a meron crystal at filling factor 
va = 0.8 in Landau level n = 0. The crystal has four merons 
per unit cell. In each unit cell, two merons with the same 
vorticity have opposite phases as explained in the text. Con- 
tours (ranging from —0.5 to 0.5) indicate the z component 
of the pseudospin with dark regions corresponding to positive 
values. 
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FIG. 2: Pseudospin texture in a meron pair crystal at filling 
factor vo — 0.8 in Landau level n = 0. The lattice is triangular 
and there are four merons per unit cell. Merons are bound 
in pairs with same value of P z and vorticities but opposite 
phases at each lattice site. Contours (ranging from —0.5 to 
0.5) indicate the z component of the pseudospin with dark 
regions corresponding to positive values. 
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FIG. 3: Hartree-Fock energy per electron as a function of 
filling factor for various crystal phases in Landau level n = 0. 



r— J— X—T, corresponding to the wave vectors {k x , k y ) = 
(0, 0) , (2tt /o) (1/V3, 1/3) , (2n/a ) (1/^3, 0) , (0, 0). 
The wave vector k represents the total distance, in 
reciprocal space and in units of 27r/ao along the path 
r — J — X — r from the origin T. The legend indicates 
in what response function x Pn , Pn , Xs*,s* , Xs y , Sy , or Xs s ,s z 
the collective mode has the biggest weight. This gives 
an indication of the nature of the mode. In Fig. |4j 
the phonon mode has its biggest weight in Xp n .p n an d 
XS^.Sx while the pseudospin wave mode has its weight 
in XSy,S y an d Xs z ,Sz- That, is, since we forced the 
pseudospin to be polarized along the x direction, the 
pseudospin wave mode corresponds to a precession of 
the pseudospin about the x axis. The phonon dispersion 
is typical of what is found for a Wigner crystal^. It is 
gapless, with u> ~ q 3 ^ 2 behavior at small wave vector. 
For the pseudospin wave, the dispersion is u> ~ q 2 at 
small wave vector confirming the SU(2) symmetry. At 
vo = 0.2, the bandwidth of the pseudospin mode is two 
orders of magnitude smaller than that of the phonon 
mode. While the bandwidth of the phonon mode does 
not change much as v$ increases to vq = 0.5, that of the 
pseudospin mode changes dramatically, becoming of the 
same order as that of the phonon mode at vq = 0.5. The 
pseudospin stiffness thus increases rapidly with vq- The 
hBCl dispersion has the same features as the eBCl as 
can be seen in Fig. \5\ 

For v > 0.55 in n — 0, we find that, within our 
numerical accuracy, the MC and MPC have the same 
energy and are lower in energy than the other phases 
considered. The dispersion relations, however, indicate 
that the MC is stable in the range v £ [0.55, 0.65] 



FIG. 4: Dispersion relation of the two Goldstone modes of the 
eBCl state at i/q = 0.2 in Landau level n = 0. The dispersion 
is plotted along the irreducible Brillouin zone of the triangular 
lattice. The left (right) y axis gives the phonon (pseudospin) 
frequency. The phonon mode has its biggest weight in Xpn-Pn 
and XSx.Sx while the pseudospin wave mode has its biggest 
weight in xs y ,s y and xs z ,s, 
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FIG. 5: Dispersion relation of the two Goldstone modes of the 
hBCl state at v\ — 0.55 in Landau level n = 1. The dispersion 
is plotted along the irreducible Brillouin zone of the triangular 
lattice. The phonon mode has its biggest weight in X Pn , Pn and 
XSx.Sx while the pseudospin wave mode has its biggest weight 
in Xs y ,s y and xs z ,s z 
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while the MPC is unstable in that range and vice versa 
for v 6 [0.65,0.90] so that there is a phase transition 
between these two states. We show in Figs. |6] and [7] 
the dispersion relations for these two states. For the 
MC phase, the dispersion is plotted along the path F — 
M — X — r, corresponding to the wave vectors (k x , k y ) — 
(0, 0) , {2w /a) (1/2, 1/2) , (27r/a ) (1/2, 0) , (0, 0) since 
the unit cell is that of a square lattice. The dispersion 
in both cases show the usual gapless phonon mode with 
lu ~ g 3 / 2 behavior at small wave vector which appears as 
a pole of Xp n ,p n an d 3 other linearly dispersing Goldstone 
modes. Some of the pseudospin modes are degenerate 
along sections of the contour of the irreducible Brillouin 
zone. The degeneracy of some of these modes is only 
lifted along r — X in the MC phase in [71 For wave vector 
k in an arbitrary direction, however, the 3 pseudospin 
modes are non degenerate. 

The energy of the two meron lattices are invariant un- 
der a rotation of the pseudospin texture around the x, y 
or z-axis since there is no equivalent of the Zeeman cou- 
pling in the graphcnc 2DEG. This implies there are three 
independent ways to rotate the spins of the state without 
any cost in energy, leading to the three Goldstone modes 
found in the GRPA. Animations of these three modes 
support this interpretation^. 




0.000 



FIG. 7: Dispersion relation of the Goldstone modes of the 
MC state at uq = 0.625 in Landau level n = 0. The dispersion 
is plotted along the irreducible Brillouin zone of the square 
lattice. The legend indicates in what response function each 
collective mode has its biggest weight. 
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FIG. 6: Dispersion relation of the Goldstone modes of the 
MPC state at vq = 0.8 in Landau level n — 0. The dispersion 
is plotted along the irreducible Brillouin zone of the triangular 
lattice. The legend indicates in what response function each 
collective mode has its biggest weight. 

We show in Fig. [H] the energies of different phases of 
the 2DEG in graphene for Landau level n = 1. We find 
the following sequence: eBCl for V\ S [0.1,0.5], hBCl 
for v Y € [0.5, 0.75] , MC for v x S [0.75, 0.80] and MPC for 
vi e [0.75,0.90]. The MC and MPC phases have almost 



the same energy within our numerical accuracy so that 
these two phases are represented by the line MC-MPC in 
[5] In comparison with the case n = 0, we see that the 
filling factor range for which a pseudospin texture exists 
for n — 1 has decreased relative to n = 0. In contrast 
with what happens in a conventional 2DEG, however, 
there is a possibility for such textures in Landau level n = 
1. This is due to the fact that the effective interactions 
in the two systems are now different in view of Eq. (|2ip . 

The dispersion of the gapless modes in Landau levels 
n = 1 and n = 2 are similar to what is seen in n = 
as can be seen from Fig. [3] Bubble crystals with more 
than one electron per site have additional gapped modes 
related to internal excitations of the bubbles; we do not 
focus on these modes in this work>£. 

Fig. [5] shows the phase diagram of the 2DEG in 
graphene for Landau level n — 2. We find here a se- 
quence of transitions involving electron and hole bubble 
crystals with one or two electrons per bubble. The two 
meron phases MC and MPC have higher energy than the 
other phases considered so that there are no meron crys- 
tals of these types in n — 2 and most probably in higher 
Landau levels as well. We emphasized however that our 
calculation is restricted to the range of partial filling fac- 
tor v — 1 > 0.1. Because the pseudospin textured states 
are pushed closer to v = 1 as n increases, we cannot rule 
the existence of meron crystals for n > 2. Indeed, Yang, 
Das Sarma and MacDonald have shown that skyrmions 
are the lowest-energy charged excitations in graphene for 
Landau levels up to n = 3^i. Meron crystals could thus 
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also be present in Landau levels n = 2,3. 

For n — 2 and near filling factor — 0.5, none of the 
phases that we considered in our analysis are stable so 
that the ground state must be of another crystal type, 
most probably the stripe phase if we compare with the 
situation in semiconductor 2DEG. This is also the con- 
clusion of Hartree-Fock calculations in Refs. dandi. A 
stripe state probably occurs near v n = 0.5 for Landau 
level n > 2 in the HFA. 



V. SIGNATURES OF WIGNER AND MERON 
CRYSTALS IN PINNING BEHAVIOR 

Generally, collective modes of 2DEG's are detected 
by inelastic light scattering or via microwave absorp- 
tion. The latter experiments are most sensitive to the 
long wavelength, low frequency behavior of the collec- 
tive modes. A disorder potential pins the Wigner crys- 
tal, in the sense that the phonon mode becomes gapped 
at a pinning frequency uj p that is dependent on the 
strength of the potential and on that of the magnetic 
field. The behavior of the pinning frequency with mag- 
netic field depends critically on the interplay between 
different length scales: the size of the electron wavefunc- 
tion on each lattice site, the lattice periodicity, and the 
magnetic lengt h 20 ' 21 . The pinning frequency in the lon- 
gitudinal conductivity may be inferred from results such 
as those found here using the replica trick^ 2 .. 

In the graphene case, both the pscudospin and phonon 
mode involve charge fluctuations so we may speculate 
that the pseudospin mode will also be pinned in the pres- 
ence of disorder, in the sense of opening a gap in their 



spectrum, because the latter generically breaks pseu- 
dospin symmetry. Disorder should also pin the four Gold- 
stone modes of the meron crystals: the phonon mode and 
the three pseudospin wave modes. If such pinned modes 
are separately observable, they could provide a unique 
signature of the formation of a Wigner or meron crys- 
tals in graphene. This conclusion should be contrasted 
with the case of a semiconductor 2DEG. There, a Wigner 
crystal has only one gapless (phonon) mode at finite Zee- 
man coupling. A skyrmion (finite Zeeman coupling) or 
meron (zero Zeeman coupling) crystal has one phonon 
mode which is gapped and one (Wigner) or three (meron) 
gapless spin wave modes that remain gapless in the pres- 
ence of disorder. 

For the pinning modes to be visible in microwave 
absorption experiments^, they must show up in 
3? [a xx (k — > 0, <J)\ where a is the conductivity tensor. 
Equivalcntly, they must appear as poles of the current - 



Xi i (k 



0,w) 



with i = 



current response functions 3 

x, y. A calculation of the conductivity tensor in the pres- 
ence of disorder is difficult for the crystal states. For a 
simple Wigner crystal (one phonon mode only), it can 
be done by mapping the system to an effective harmonic 
mode l 20 ' 22 and using the replica trick. So far there has 
been no generalization of this method to crystals with 
an additional layer or valley degree of freedom. To test 
our speculation that the phonon and valley pseudospin 
modes are visible in microwave absorption, we "simulate" 
a disorder potential by adding a periodic external poten- 
tial. Note that our HFA method forces us to choose this 
periodicity to be the same as that of the crystal consid- 
ered. 
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In formulating the relevant response function for con- 
ductivity, the current operator is built up from opera- 
tors that excite electrons between Landau levels. This 
presents problems when the Hilbert space is restricted 
to one Landau level as in our calculation. In principle 
one needs to retain other Landau levels in order to ob- 
tain a non- vanishing result, significantly complicating the 
calculations. To circumvent this, we wish to find an ap- 
propriate projection of the current operator into a sin- 
gle Landau level. To do this, we generalize a method 
first introduced by Girvin, MacDonald and Platzman 
in Ref. This procedures captures the drift current 
j (r) = ^p(r) VJ7 (r) x z of the electrons in the poten- 
tial U (r) and gives a current that satisfies the continuity 
equation. Improvements upon this procedure are possi- 
ble but lead to very complicated expressions 2 — . We first 
write a second quantized hamiltonian including the (to- 
tal) density p (k) and valley-pseudospin operator P (k) 



H = 



1 f e 



4Uf 



-No.. 



£T(q)p(-q)p(q) (48) 
q 

£X(q)P(-q)-P(q), 



where 



T(q)=2tf(q)-X(q). 



(49) 



With the hamiltonian of Eq. (I48[) , we obtain the equation 
of motion of the density operator 



^ = [p(k) )ff ], 



(50) 



which we linearize by writing p(k) — > (p(k)) + Sp (k) 
where the average is evaluated in the HFA. Keeping terms 
up to linear order in Sp (k) ,we find 



n 



(k) 



dt 



iL)£[T(G)-T(k-G)] 



(51) 



G 



x sin (k x Gf/2) (p (G)) Sp (k - G) 

- 4 (5)|> (G) - X(k - G)] 

x sin (k x G£ 2 /2) (P (G)) • 5P (k - G) . 

This is the equation of motion of the density in the 
GRPA. 

To find an expression for the current operator valid at 
small k, we make the approximation 



sin (k x G£ 2 /2) m k- (G x z£ 2 /2) 
and use the continuity equation 
dSp (k) i 



dt 



-k-j(k) 



(52) 



(53) 



We get in this way 



: £(Gx») [T (k + G)-T(G)] 



G 



(p(-G))Sp(k + G) 



(54) 



+2 Mr E( Gx *0 (k + G) - X (G)] 



tlK 

x (P(-G))-<5P(k + G). 

This expression shows that the current can have contri- 
butions from both density and pseudospin fluctuations. 
In terms of the original pij (q) operators, we can write 
the current expression as 



j (k) = -i— Y, E ( k < G ) «W (k + G) , (55) 
where 



G a,b 



F a , fc (k,G) (56) 
= (G£ xz)[H(k + G)-H(G)](p(-G))S a , b 
- (G£ xz)[I(k + G)-I (G)] (p M (-G)) . 

With Eq. (p)5"|) for the current, we can easily write the 
current-current Matsubara Green's function tensor 



X J < J (k,T)--A^(Ti(k,T)j(-k,0)> 



(57) 



so that the retarded current-current response function is 
finally given by 



(58) 



„b,c,dG,G' 



x x a n Ac - d (k + G, k + G', lu) ¥ c . d (-k, -G') . 



If we apply an external potential, the hamiltonian H 
H + H ext with 



H ext = nJ — ) Y E W « (q) 



(59) 



a q 



We allow the potential W a (q) to be different for the 
K and K' valleys. With this potential, the function 
F aj f, (k, G) in Eq. ([58]) must be replaced by 

F id (k, G) (60) 
= - (G x zi) W K (-G) S id 5i, K 

— (G x zl) Wk> (-G) 

+ (G£ xS)[ff(k + G)-ff (G)] (p (-G)) S itj 

- (G£ xz)[I(k + G)-I (G)] (pj,i (-G)) . 

In pseudospin language, this means that the current in 
Eq. (f5"4"| becomes j (k) — > j (k) + }w (k) , where 



jw (k) 



4 \M 



2\hn 



E (G x zl) W+ (-G) Sp (k + GJB1) 

G 

E (G x z£) W- (-G) 5P Z (k + G) , 
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where W ± (-G) = W K (-G) ± Wjc (-G) . Note that 
for consistency, we also include the external potential in 
the calculation of the (pi j (G)) s as well as in that of 
^(k + G.k + G'^)'. 

In the absence of any external potential, we find that, 
for the Wigner or meron crystal, only the phonon mode 
(and some higher energy modes) appears as a pole of 
3 [Xxx for any wave vector k. The phase modes 

are conspicuously absent of the current response. We 
explain this by the fact that the phase modes are trans- 
verse modes: the motion of a pseudospin is perpendicu- 
lar to the local value of that pseudospin so that the term 
(P(-G')) ■ <5P(k + G') - in Eq. {54]). Also, these 
modes have no weight in the density-density response 
function. They cannot contribute to the local current. 
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FIG. 10: Dispersion relation of the phonon and pseudospin 
modes of the Wigner crystal at vq — 0.4 in Landau level 
n = in an external potential Wk = —0.005 and Wk' = 0. 
The dispersion is plotted along the irreducible Brillouin zone 
of the triangular lattice (see Fig. 5). Both the phonon and 
pseudospin wave mode are gapped by the external field. 

This conclusion is unchanged, for the Wigner crystal 
state, if we apply an external potential with Wk (q) = 
Wk> (q) • The phonon mode is gapped by that poten- 
tial but the pseudospin mode dispersion remains gap- 
less because the pseudospin symmetry has not been bro- 
ken. To induce a gap in the pseudospin mode, we must 
allow Wk (q) ^ Wk 1 (q) so that the external poten- 
tial can couple to the z component of the pseudospin: 
Pn,z (q) = (q) - Pn'~ (q))/2- We show in Fig. 10 

the dispersion relation of the collective modes of the 
Wigner crystal for v = 0.4 and with an external potential 
that is different in the two valleys i.e. Wk — —0.005 and 
Wk> = 0. Both the phase and phonon modes are now 
gapped as expected but Fig. 11 shows that, once again, 



only the phonon appears as pole of the current. This fact 
can be readily understood: the external potential is much 
weaker than the exchange energy that forces the parallel 
alignment of the pseudospins. It follows that when the 
external potential is applied, the pseudospins all align 
along the z axis even though the external potential is 
modulated in space. The pseudospin mode is a trans- 
verse mode and so both terms (P (— G')) • SP (k + G') 
and 5P Z (k + G) in the definition of the current in Eq. 
(f5"4"f are zero. For the phonon mode in the Wigner crys- 
tal state, the peak in the current response stays finite 
with decreasing wave vector so that the phonon mode is 
visible at k = and can contribute to the absorption. 
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FIG. 11: Imaginary part of the current response function. 
From left to right: k x = k y = 0.001, 0.01, 0.05, 0.1, 0.3 in units 
of 2n/ao- These peaks come from the phonon mode. There is 
no peak coming from the pseudospin mode. 

For the meron crystal, the situation is more compli- 
cated. An external potential W- (— G) acts as a pseudo- 
magnetic field. If it is uniform in space, it acts effectively 
as a pseudospin Zeeman coupling and produces a transi- 
tion from a meron to a bimeron crystal with two Gold- 
stone modes (the phonon mode and a gapless pseudospin 
mode related to the U(l) symmetry of the hamiltonian) 
and two gapped pseudospin modest. The pseudospin 
modes being transverse and the second term in }w (k) 
being zero because of the uniformity of the potential, 
there should again be no contribution of these modes to 
the conductivity. This is the case for two of the pseu- 
dospin modes, but the gapless pseudospin mode has a 
weight in the density-density response function at finite 
wave vector k and does appear in the current response 
along with the phonon peak. Both peaks go to zero with 
decreasing wave vector however. 

For the meron crystal, the second term in jw (k) 
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should be finite if W_ (— G) is non-uniform, in contrast 
with the Wigner crystal case, because some of the pseu- 
dospin modes involve a fluctuation in the z component of 
the pseudospin. Unfortunately, we find the meron crystal 
to be very sensitive to an inhomogcneous external poten- 
tial, as indicated by instabilities in the collective mode 
spectrum at small values of k. This is likely related to 
the extreme closeness in energy of the meron crystal and 
meron pair crystal states, so that the external potential 
may lead to a different and possibly more complicated 
textured state. Thus it is not possible to compute the 
response functions at small wave vector k for these tex- 
tured states without better knowledge of their ground- 
state structure in the presence of a pinning potential. 
We note that at finite wavevector, where the dispersion 
is well behaved, the pseudospin modes do appear in the 
current response as expected. 

In spite of the difficulty demonstrating the presence of 
a signature of the pinned phase modes in the dynamical 
conductivity at small wavevector, we believe at least a 
small response will in fact generically always be present. 
Beyond the density fluctuations in the pseudospin modes 
due to spin-charge coupling, there is a further density 
response due to the fact that the A and B sites of the 
lattice are at difference positions in real space in a unit 
cell. In our calculations, A and B were treated as two 
orthogonal "spin" states of an electron, but their slightly 
different locations in real space were not included in the 
model. If this were included, we expect that an oscilla- 
tion of the pseudospin that changes the relative weight 
of an electron on the A or B sublattice (or, equivalently, 
on the K and K' valleys) will translate into a change in 
the position of that electron or into a dipole fluctuation. 
Thus, if this distinction were properly included in our 
model, we would expect that the phase mode would ap- 
pear as a pole of the current response just as the phonon 
mode does, albeit weakly, since the symmetry breaking 
is small. 

In closing this section, we remark that, in Fig. [H the 
effective stiffness for the pseudospin mode is two orders 



of magnitude smaller than that of the Wigner crystal for 
filling factor v < 0.2. In the presence of disorder, we 
might then expect two pinning modes of very different 
frequencies for a Wigner crystal and it may be impossi- 
ble to detect the two modes simultaneously in an actual 
experiment. The meron crystal dispersion does not suf- 
fer from this problem since all four modes appear to have 
similar band widths. 



VI. CONCLUSION 

We have shown in this work that the 2DEG in graphcnc 
can support Wigner crystals and meron crystals with 
valley-pseudospin textures. Our numerical analysis was 
restricted to filling factor v n £ [0.1,0.9] in each Landau 
level n and we concluded that, in this range, meron crys- 
tals are present in Landau levels n — and n — 1 only. 
We have computed the dispersion relation of the collec- 
tive excitations of these two crystal states and showed 
that the Wigner crystal has one extra Goldstone mode 
with a quadratic dispersion at small wave vector in ad- 
dition to the phonon mode. Meron crystals have 3 ex- 
tra Goldstone modes in addition to the phonon mode. 
These extra Goldstone modes are valley-pseudospin fluc- 
tuations. In graphene, all these modes involve density 
fluctuations, and we speculated that these last modes 
could be visible as pinning modes in microwave absorp- 
tion spectrum in a real disordered system. 



Acknowledgments 

This work was supported by a research grant from the 
Natural Sciences and Engineering Research Council of 
Canada (NSERC) for R. Cote and an NSF Grant No. 
DMR-0704033 for H. A. Fertig. Computer time was pro- 
vided by the Reseau Quebecois de Calcul Haute Perfor- 
mance (RQCHP). 



1 P. K. Lam and S. M. Girvin,Phys. Rev. B 30, 473 (1984); 
D. Levesque, J. J. Weis, and A. H. MacDonald, Phys. Rev. 
B 30, 1056 (1984); K. Esfarjani and S. T. Chui, Phys. Rev. 
B 42, 10758 (1990); K. Yang, F. D. M. Haldane, and E. H. 
Rezayi, Phys. Rev. B 64, 081301 (R) (2001); X. Zhu and 
S. G. Louie, Phys. Rev. B 52, 5863 (1995). 

2 For recent reviews, see Physics of the electron solid, edited 
by S. T. Chui (International Press, Boston (1994) and H. 
Fertig and H. Shayegan in Perspectives in Quantum Hall 
Effects, edited by S. Das Sarma and A. Pinczuk (Wiley, 
New York, 1997), Chaps. 5 and 9 respectively. 

3 P. D. Ye, L. W. Engel, D. C. Tsui, R. M. Lewis, L. N. 
Pfeiffer, and K. West, Phys. Rev. Lett. 89, 176802 (2002); 
Yong P. Chen, G. Sambandamurthy, Z. H. Wang, R. M. 
Lewis, L. W. Engel, D. C. Tsui, P. D. Ye, L. N. Pfeiffer, 
and K. W. West, Nat. Phys. 2, 452 (2006); Y. P. Chen, R. 



M. Lewis, L. W. Engel, D. C. Tsui, P. D. Ye, L. N. Pfeiffer, 
and K. W. West, Phys. Rev. Lett. 91, 016801 (2003); R. M. 
Lewis, Yong Chen, L. W. Engel, D. C. Tsui, P. D. Ye, L. 
N. Pfeiffer, and K. W. West, Physica E22, 104 (2004); R. 
M. Lewis, P. D. Ye, L. W. Engel, D. C. Tsui, L. N. Pfeiffer, 
and K. W. West, Phys. Rev. Lett. 89, 136804 (2002); R. 
M. Lewis, Y. Chen, L. W. Engel, D. C. Tsui, P. D. Ye, L. 
N. Pfeiffer, and K. W. West Phys. Rev. Lett. 93, 176808 
(2004); R. M. Lewis, Yong Chen, L. W. Engel, P. D. Ye, 
D. C. Tsui, L. N. Pfeiffer, and K. W. West, Physica E22, 
119 (2004); Yong P. Chen, R. M. Lewis, L. W. Engel, D. 
C. Tsui, P. D. Ye, Z. H. Wang, L. N. Pfeiffer, and K. W. 
West, Phys. Rev. Lett. 93, 206805, 2004. 
4 A. A. Koulakov, M. M. Fogler and B. I. Shklovskii, Phys. 
Rev. Lett. 76, 499 (1996); M. M. Fogler, A. A. Koulakov, 
and B. I. Shklovskii, Phys. Rev. B 54, 1853 (1996); R. 



14 



Moessner and J.T. Chalker, Phys. Rev. B 54, 5006 (1996); 
M. M. Fogler and A. A. Koulakov, Phys. Rev. B 55, 9326 
(1997). For a review of the bubble and stripe phases in 
higher Landau levels, see M. Fogler in High Magnetic 
Fields: Applications in Condensed Matter Physics and 
Spectroscopy, ed. by C. Berthier, L.-P. Levy, G. Martinez 
(Springer- Verlag, Berlin), 99 (2002). 

R. Cote, C. B. Doiron, J. Bourassa, and H. A. Fertig, Phys. 
Rev. B 68, 155327 (2003). 

C.-H. Zhang and Yogesh N. Joglekar, Phys. Rev. B 75, 
245414 (2007). 

C.-H. Zhang and Y. N. Joglekar, Phys. Rev. B 77, 205426 
(2008). 

M. I Katsnelson, K. S. Novoselov, and A. K. Geim, Nature 
Physics 2, 620 (2006). 

J. Wang, H.A. Fertig, A. P. Iyengar, and L. Brey, unpub- 
lished. 

M. O. Goerbig, R. Mosessner, and B. Dougot, Phys. Rev. 
B 74, 161407(R) (2006). 

K. Yang, S. Das Sarma, and A. H. MacDonald, Phys. Rev. 
B 74, 075423 (2006). 

L. Brey, H.A. Fertig, R. Cote, and A.H. MacDonald, Phys- 
ica Scripta, T 66, 154 (1996). 

A. H. Castro Neto, F. Gui nea, N. M. R. P eres, K. S. 
Novoselov, and A. K. Geim, larXiv:0709.TT63l (To be pub- 
lished in Review Modern Physics (2008)). 



1 J. Bourassa, B. Roostaei, R. Cote, H. A. Fertig, and K. 
Mullen, Phys. Rev. B 74, 195320 (2006). 

15 W. Luo and R. Cote, unpublished. 

16 R. Cote, D. B. Boisvert, J. Bourassa, M. Boissonneault, 
and H. A. Fertig, Phys. Rev. B 76, 125320 (2007). 

17 R. Cote and A. H. MacDonald, Phys. Rev. B 44, 8759 
(1991); Phys. Rev. Lett. 65, 2662 (1990). 

18 A. P. Iyengar, J. Wang, H.A. Fertig, and L. Brey, Phys. 
Rev. B 75, 125430 (2007). 

19 Y.V. Nazarov and A.V. Khaetskii, Phys. Rev. Lett. 80, 
576 (1998). 

20 R. Chitra, T. Giamarchi, and P. Le Doussal, Phys. Rev. 
B65, 035312 (2001); ibid., Phys. Rev. Lett. 80, 3827 
(1998). 

21 H. A. Fertig, Phys. Rev. 59, 2120 (1999). 

22 R. Cote, Mei-Rong Li, A. Faribault, and H. A. Fertig, 
Phys. Rev. B 72, 115344 (2005). 

23 S. M. Girvin, A. H. MacDonald, and P. M. Platzman, Phys. 
Rev. B 33, 2481 (1986). 

24 J. Martinez and M. Stone, International Journal of Mod- 
ern Physics B 7, 4389-4401 (1993); R. Rajaraman, Inter- 
national Journal of Modern Physics B 8, 777-788 (1994); 
R. Rajaraman and S. L. Sondhi, Modern Physics Letters 
B 8, 1065-1073 (1994). 



